Anomalous avoided level crossings in a Cooper-pair box spectrum 
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Abstract 

We have observed a few distinct anomalous avoided level crossings and voltage dependent tran- 
sitions in the excited state spectrum of an Al/A10x/Al Cooper-pair box (CPB). The device was 
measured at 40 mK in the 15 - 50 GHz frequency range. We find that a given level crosses the CPB 
spectrum at two different gate voltages; the frequency and splitting size of the two crossings differ 
and the splitting size depends on the Josephson energy of the CPB. We show that this behavior is 
not only consistent with the CPB being coupled to discrete charged "two-level" quantum systems 
which move atomic distances in the CPB junctions but that the spectra provide new information 
about the fluctuators, which is not available from phase qubit spectra of anomalous avoided levels. 
In particular by fitting a model Hamiltonian to our data, we extract microscopic parameters for 
each fluctuator, including well asymmetry, tunneling amplitude, and the minimum hopping dis- 
tance for each fluctuator. The tunneling rates range from less than 3.5 to 13 GHz, which represent 
values between 5% and 150% of the well asymmetry, and the dipole moments give a minimum 
hopping distance of 0.3 to 0.8 A. We have also found that these discrete two-level systems have a 
pronounced effect on the relaxation time (Ti) of the quantum states of the CPB and hence can be 
a source of dissipation for superconducting quantum bits. 

PACS numbers: 85.25. Cp and 85.35. Gv 
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Two-level fluctuators (TLF) of one kind or another have long been recognized as an un- 
derlying cause of resistance fluctuations in metals,- 1 ^ critical current fluctuations in Joseph- 
son junctions,— excess flux noise in SQUIDs,-^ and telegraph noise and excess charge 
noise in Coulomb blockade devices.— For charge and critical current fluctuations - 
the fluctuators appear to be moving charges or rotating electric dipoles in the insulating 
junction barrier or nearby dielectric layers 3 - but the precise microscopic identification of 
the fluctuators is not settled. 

There are several reasons why it has been difficult to be certain of the precise microscopic 
agents causing charge or critical current fluctuations. First, the effects from individual 
fluctuators are typically very small, especially for the cryogenic or milli-Kelvin temperatures 
of interest here, and this typically makes the resulting experiments quite challenging. Second, 
a variety of materials and processes have been used to build devices, and the presence of 
a fluctuator may well depend on both the materials and the fabrication technique. Third, 
while measurements of telegraph noise can reveal microscopic information about individual 
fluctuators, in many cases it is possible only to distinguish the largest one or two fluctuators 
in a background noise of smaller fluctuations. Compared to the number of atoms in even 
the smallest device, such observable discrete fluctuators are extremely rare and thus may 
not be representative of a typical atomic scale defect in the device. Finally, much of the 
experimental data obtained on fluctuators consists of relatively smooth 1/f noise power 
spectra. Such noise spectra arise from a distribution of many fluctuators^ 1 ^ and it is often not 
possible to resolve individual fluctuators. From smooth spectra, it is difficult to determine 
uniquely such basic microscopic parameters as the hopping distance or the absolute number 
of fluctuators in a given energy range. 

Research in quantum computing based on superconducting devices has led to increased 
interest in understanding two-level fluctuators. In qubit research, fluctuators can be a se- 
rious problem because they can lead to decoherence, dissipation, inhomogenous broadening 
and a decrease in measurement fidelity— ii£»iZii2ai2i2Jl With this new interest have also come 
new approaches to the problem. In particular, microwave spectroscopy of the various su- 
perconducting qubits has revealed the presence of small un-intended avoided crossings in 
the transition spectrum due to coupling between the device and individual two-level fluctu- 
ators.— ii£iiSiiSi2L2i Avoided level crossings were initially seen in the phase qubit^ and sub- 
sequently in the flux qubit 18 ' 20 , quantronium^, and the transmon.— While it was initially 
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thought that the anomalous avoided level crossings were due to critical current fluctua- 
tions,— a later comparison of a hopping distance (extracted from an analysis of an ensemble 
of avoided level crossings) to atomic distances suggested that charge fluctuators were re- 
sponsible.— 

In this article, we report observations of anomalous avoided level crossings in the tran- 
sition spectrum of a Cooper-pair box (CPB), a superconducting quantum device that is 
sensitive to charge. We find that each avoided level crossing is due to a charged two level 
system in which the charge moves atomic distances in the tunnel barrier of the device's 
Josephson junction. The spectra contain a striking feature not seen in phase qubits - the 
two level system's levels are voltage dependent - and by modeling the full spectrum, we can 
resolve some of the key microscopic parameters of each fluctuator. We also find a strong 
correlation between the locations of the avoided level crossings and decreases in the lifetime 
(Ti) of the excited state of the qubit. 

For our measurement, we use a CPB consisting of a small superconducting island con- 
nected to superconducting leads by two ultra-small (nominal area 120 nm by 120 nm) 
Josephson junctions (see Fig. CQ. By applying voltage V g to a gate electrode, which is 
capacitively coupled to the CPB island with capacitance C g , we can control the number n of 
excess Cooper pairs on the island. Neglecting the quasiparticle states,— 1 ^ the Hamiltonian 
describing the CPB is given by2^2£ 

Hcpb = E c £(2n - n g f \n) (n\ - ^ + 1) (n\ + \n)(n + l\) (1) 

where n g = —C g V g /e is the reduced gate voltage, e = —1.609 x 10~ 19 C is negative, E c = 
e 2 /2Cs is the charging energy, Cs is the total capacitance of the island, and Ej is the 
Josephson energy. At the CPB charge degeneracy point, which occurs when n g is an odd 
integer, the energy difference between the first excited state and ground state is a minimum 
and approximately Ej. By applying a magnetic flux (0) through the superconducting loop, 
Ej can be effectively changed by the quantity Ej = Ej Iax cos(7T0/</> o ) where Ej iax = hl c /2e, 
I c is the critical current of the two junctions, and O = hc/2e is the magnetic flux quantum. 
Similarly, the first and second excited state have a minimum splitting of approximately 
Ej/1QE C and this minimum occurs at even values of n g . 

To measure the excess charge on the CPB, we used a superconducting Coulomb-blockade 
electrometer that is capacitively coupled to the CPB [see Fig. Ufa)].— The electrometer 
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and the CPB were formed together on a quartz substrate in the same process; we used 
electron beam lithography and standard double-angle evaporation of Al with an oxidation 
step to form an AlOx tunnel barrier between the two Al layers. 27 To improve the bandwidth 
of the measurement, the electrometer was connected to an on-chip LC tank circuit and 
operated in an rf-SET mode;— the reflectance was measured at the resonance frequency of the 
tank circuit (640 MHz) while the electrometer was dc biased at the Josephson quasiparticle 
resonance and with the gate of the electrometer tuned to maximize the sensitivity to charge 
changes. ~ 

The sample was mounted in a Cu box that was bolted to the mixing chamber of an Oxford 
Instruments model 100 dilution refrigerator with a base temperature of 40 mK. Figure [2] 
shows a plot of the measured change in the rf reflectance of the rf-SET as a function of n g 
and frequency of a microwave perturbation applied to the gate of the CPB. For this plot, 
the microwave frequency was adjusted between 24 and 50 GHz in steps of 30 MHz. When 
the applied frequency is resonant with a level splitting of the CPB, the CPB state changes, 
causing a change in the charge on the island, which causes a change in the reflectance of 
the rf-SET. The white parabolic like band between < n g < 1 in Fig. [2] corresponds to 
a measured change of ~ le on the island of the CPB while the black parabola between 
1 < n g < 2 corresponds to a measured change of ~ — le on the island of the CPB. 

Close examination of the spectrum reveals a few small avoided crossings which are not 
predicted by Eq. [TJ anomalous avoided crossings imply that additional degrees of freedom 
are coupled to the CPB. Figure [3] shows a detailed view of the spectrum near two prominent 
avoided crossings due to a single fluctuator for four different values of Ej. In Fig. [3](a), one 
crossing occurs at / = 34.3 GHz and n g = —0.43. A second crossing occurs at / = 33.5 
GHz and n g = 0.48. We note in particular the presence near each crossing of short features 
in the spectrum that appear to point toward the other crossing (see arrows in Fig. 3a which 
indicate small projections in the spectrum), suggesting the two crossings are related. When 
n g was swept over multiple periods, we found that this spectrum was periodic with period 
2. Figures [3](a) through [3](d) show that the spectrum changes as Ej/ks is decreased from 
1.0 K to 0.1 K. In particular as Ej is reduced, the avoided crossing splitting decreases from 
a size of 150 MHz to a size that is too small for us to determine. The other anomalous 
avoided crossings we observed showed similar behavior. 

The dependence of the crossings on gate voltage and Josephson energy are consistent with 
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the CPB being coupled to a charged fluctuator that is moving in one of the Al/A10x/Al 
ultra-small junction barriers. In particular, because the avoided level crossings have a peri- 
odicity of 2 in reduced gate voltage, this suggests that the charge fluctuator resides in one of 
the two tunnel junctions that form the CPB; the electrostatic potential of the CPB island 
is given by Vi = e/Cs(2n — n g ) which has a periodicity of 2 in n g for the different states of 
the CPB. 

In general, the Hamiltonian for a charged fluctuator would depend on the local environ- 
ment including other charged defects in the system. To simplify the problem, we assume 
that the fluctuator takes the form of a point charge moving between two potential wells. In 
this approximation the Hamiltonian for the TLF is 

Htlf = £ a \x a ) (x a \ + e b \x b ) (x b \ + T ab (\x a ) (x b \ + \x b ) (x a \) (2) 

where e a is the energy of the fluctuator at position x a , e b is the energy of the fluctuator 
at position x b , and T ab is the tunneling matrix element between the two sites. For an 
isolated fluctuator, the difference in energy between the two states of the TLF is given by 



(e b - e a y + 4T a 2 6 . 

Assuming the fluctuator is a point charge that resides in the tunnel junction, with charge 

Qtlf-i the change in the induced polarization charge on the island of the CPB when the 

fluctuator tunnels from position x a to x b is given by AQ pi = Qtlf^ — x a ) cos(rj)/d where r\ 

is the angle of displacement relative to the electric field in the junction and d is the thickness 

of the tunnel junction.— Since the charge on the island is constant when n is constant, a 

corresponding change in the polarization charge (—AQ pi ) must exist over the rest of the 

island. This in turn induces a change in the induced surface charge on the gate electrode 

given by AQ g = AQ^CgjC-^- Finally, we find that the change in the electrostatic potential 

of the island due to motion of the point charge is given by^ 

AVi = Qtlf (x b - Xg) cos(ry) 
Cs d 

Note that the total electrostatic potential of the island depends on the number of excess 
Cooper-pairs on the island, the displacement of the point charge and the reduced gate 
voltage. Accounting for the electrostatic charging energy and the work done by gate voltage 
source when the point charge moves, the following coupling Hamiltonian is found^ 

H CPB _ TLF = 2E C 2^ 2J2n-n 9 ) \n) \x) {x\ (n| . (4) 

x=x a ,xi, n C & 



Combining Eqs. [TJ [2j and H] we find the total Hamiltonian for a CPB coupled to a single 
fluctuator. From this, we can find the energy levels and the transition frequencies from the 
ground state to the excited states of the system and fit the measured spectrum. We note that 
in this model, an avoided crossing occurs when the first excited state of the TLF is resonant 
with the first excited state of the CPB. The change in the induced polarization charge on 
the CPB island due to the fluctuator being excited ultimately gives rise to the avoided level 
crossings at two different frequencies and two different reduced gate voltages (i.e. it breaks 
the symmetry of the CPB). The splitting size in our model depends on Ej; coupling the two 
excited states together is a second order process that involves both the tunneling of a Cooper- 
pair and the tunneling of the TLF. The dashed red curve in Fig. [3] shows the predicted 
spectrum for the splitting we found near 34 GHz. Here we used fit parameters (e& — e a )/kB 
= 1.427 K, T a b/kB = 0.39 K, and Qtlf^x cos(r]) Jed = —0.074 where Sx = (x\, — x a ) is the 
maximum displacement of the fluctuator. We note that extracting these three parameters 
requires us to measure the splitting size and frequency of the avoided crossings at two 
different gate voltages and that these fit parameters can be varied by approximately 10% to 
maintain a good fit. 

We found another prominent avoided crossing near 43 GHz and smaller avoided crossings 
near 20 and 39 GHz. Table [J summarizes the best fit results for all of the observed avoided 
crossings. The avoided crossing near 39 GHz had a splitting size too small to resolve, which 
places an upper bound on T a b of the TLF. Assuming a TLF charge of Qtlf = e and a tunnel 
barrier thickness of d = 1 nm, we extract minimum hopping distances for the fluctuators 
that range from 0.32 A to 0.83 A. 

The model also predicts state transitions of the TLF at gate voltages away from the 
observed avoided crossings. Although these were not visible in Fig. [31 we found that weak 
transitions could be observed far from the avoided crossing when the microwave excitation 
power was increased by approximately a factor of ten. Figure H](a) shows the measured 
excess charge spectrum between an applied frequency of 33 and 34.6 GHz; note the very 
faint transition due to the TLF between n g = -1.25 and n g = -0.75. In this region the 
measured excess charge on the CPB is approximately -0.03 e, about an order of magnitude 
smaller than predicted from our simple theory. Figure H](b) overlays the predicted spectrum 
using the parameters found in Fig. [3j 

Similarly, Fig. H^c) shows the measured excess charge spectrum between 41.5 to 44.5 GHz; 
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a frequency range where another prominent avoided crossing was observed. Note that in 
Fig. H](c) we observe weak voltage dependent transitions, due to the TLF, near n g = 0.8 and 
/ = 44.25 GHz; as well as between n g = —0.8 and -0.4. Again the predicted spectrum in 
Fig. SJd) follows the measured spectrum accurately. In this spectrum, we observe additional 
features like non-equilibrium quasiparticles which creates the apparent periodicity in n g of 
2 22j23 mu ltiple CPB spectra (i.e. spectra parallel to the CPB spectrum) which is interpreted 
to be due to low frequency charge ffuctuators, and horizontal bands which are due to a change 
in the gain of the Coulomb-blockade electrometer. 

We have also measured the lifetime (Ti) of the excited state of the CPB as a function of 
the transition frequency. For these measurements, the CPB was excited at its zero to one 
transition frequency, and the charge on the island was continuously monitored as a function 
of time using the rf-SET after the excitation source was turned off. The measurements were 
done at a small Josephson energy to decouple the CPB from charge perturbations in the 
system and hence increase the lifetime of the qubit to a maximum value.— » Figure [5] shows 
the measured decay rate (Ti) from 15 GHz to 45 GHz at Ej/ks = 0.1 K. Several peaks in 
the decay rate are evident; we find that when the CPB is in resonance with the TLF, the 
measured lifetime decreases from a few microseconds to 1 /xs or less. This decrease in 7\ 
near a resonance was a useful tool for finding some of the avoided crossings. This behavior 
suggests that the lifetime of these TLFs is smaller than a few microseconds and that the 
interaction of the CPB with a charged TLF is a source of dissipation for the CPB. 34 

Finally, to investigate the stability of the TLFs, the device was warmed up to room 
temperature. After "annealing" at room temperature for 14 days, the device was cooled 
again to 40 mK. We found an avoided crossing associated with one level in the 20 - 50 GHz 
range and one peak in the decay rate at the same frequency. The new crossings occurred 
around / = 23 GHz and fitting to the charge model gave: (e& — e a )/kB = 0.34 K, Tab/ks 
= 0.52 K, and Qtlf^x cos(rj) fed = 0.078. We also note that the fluctuators observed here 
were stable after annealing at 4 K for two days or after placing the device in the normal 
state by applying a 1 T magnetic field for one hour. 

In conclusion, we have observed unintended voltage-dependent transitions and avoided 
level crossings in the excited state spectrum of a CPB, consistent with a charge fluctuator 
that can tunnel between two locations separated by atomic scale distances in the tunnel 
junction. Finally, we note that the spectra allow us to extract key microscopic parameters 
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of the TLF, such as the tunneling matrix element, and test some theories of fluctuators in 
these devices.- 1 ^ 1 ^ 
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FIG. 1: (a) Schematic of experimental setup. A rf Coulomb-blockade electrometer measures the 
charge state and energy spectrum of the Cooper-pair box (CPB). The CPB and electrometer had 
a charging energy of E c /kB=0-575 K and E Ct E/kB=l-l K, respectively, (b) Scanning electron 
micrograph of the device. 

FIG. 2: Measured spectrum of CPB when the applied microwave frequency / is adjusted from 
24 GHz to 50 GHz with Ej/kg= 1.12 K. The white parabolic like band between < n g < 1 
corresponds to a measured change of ~ le on the island of the CPB while the black parabola 
between 1 < n g < 2 corresponds to a measured change of ~ — le on the island of the CPB. 

FIG. 3: Measured spectrum of CPB around 34 GHz at different Josephson energies as specified 
in the graphs. The red and blue colors represent the charge on the island corresponding to the 
excited and ground state of the system, respectively. The arrows in (a) indicate a small projection 
of the avoided level that appear to point toward one another. The red dashed line is the predicted 
spectrum using a Hamiltonian consisting of a charged two level fluctuator (see Table U for the 
microscopic parameters of the TLF) coupled to a CPB. 

FIG. 4: Measured excess charge spectrum at relatively large microwave drive amplitudes near two 
of the more prominent avoided level crossings, (a) Measured spectrum between 33 and 34.6 GHz 
with Ej/kB= 0.95 K (same avoided level crossing as in Fig. [3]). (b) Same as (a) with the predicted 
spectrum from our system Hamiltonian (red dashed curve) using the parameters in Fig. [3] and 
Table U and observed local peaks in the spectrum due to the TLF as the blue points, (c) Measured 
spectrum between 43 and 44.5 GHz at Ej/kB= 0.1 K. (d) Same as (c) with predicted spectrum 
(red dashed curve) using the parameters in Table [J and observed local peaks in the spectrum due 
to the TLF as the blue points. 

FIG. 5: (Color online) Measured decay rate (Ti) at Ej/ks = 0.1 K as a function of energy 
level separation after first cooldown (blue triangles) and after the device was warmed up to room 
temperature and cooled back down to T = 40 mK (red circles). The observed peaks in T\ before 
annealing (20, 34, 39, and 43 GHz) and after annealing (23 GHz) correlate with the location of 
observed avoided level crossings (see Table H]). 
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TABLE I: Parameters of observed TLFs. N is the TLF number, / is the frequency near where the 
avoided level crossing occurs, n g ^ + is the negative and positive reduced gate voltage location of 
each TLF and A~/+ represents the splitting size of the crossing at n g ^ + for the specified value Ej. 
(eb — e a ) is the asymmetry in the potential energy of the TLF in the two wells, Q5x cos(ry) / 'ed is the 
coupling between the TLF and the CPB and T a b is the tunneling term between the two potential 
wells. 



N 


f (GHz) 


n g 


A" (MHz) 


n+ 


A+ (MHz) 


Ej/k B (K) 


{£b -£a)/k B (K) 


Qtlf&x cos(rj)/ed 


T ab /k B (K) 


1 


20 


-0.60 


30 


0.63 


10 


0.3 


0.75 


-0.032 


0.3 


2 


34 


-0.43 


140 


0.48 


80 


1.0 


1.43 


-0.074 


0.39 


3 a 


39 


-0.34 




0.36 




1.1 


1.86 - 1.89 


-0.083 


< 0.17 


4 


43 


-0.26 


120 


0.25 


220 


1.1 


1.61 


0.083 


0.65 


5 b 


23 


-0.57 


90 


0.58 


140 


0.55 


0.34 


0.078 


0.52 



a The splitting size was too small to resolve. 
^After annealing device at T = 300 K for 14 days. 
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